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A quadrupole pattern of the out-of-plane component of the magnetic field inside a reconnection 
! region is seen as an important signature of the Hall-magnetohydro dynamic (Hall-MHD) regime of 
■ reconnection. It has been first observed in numerical simulations and just recently confirmed in 
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. the MRX (Magnetic Reconnection Experiment) [M. Yamada, H. Ji, S. Hsu, T. Carter, R. Kulsrud, 
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field is analyzed and traced to a current of electrons that flows along the lines in and out of the 



^ ' inner reconnection region to maintain charge neutrality. The role of the quadrupole magnetic field 
in the overall dynamics of the reconnection process is discussed. In addition, the bipolar poloidal 
electric field is estimated and its effect on ion motions is emphasized. 
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I. INTRODUCTION 



Ever since it was established that the classical Sweet-Parker reconnection model 2j, |3j, |4 1 
with Spitzer resistivity is too slow and that the Petschek |5J| fast-reconnect ion mechanism 
cannot be realized in resistive MHD with uniform resistivity 6 ), theoretical stud- 

ies of fast reconnection have proceeded as a competition between two schools of thought. The 
first one invokes the idea of enabling the Petschek mechanism by a strongly localized anoma- 
lous resistivity due to 



asma micro-instabilities triggered when a certain current threshold 
121 ] . The second, commonly referred to as the Hall reconnection mech- 
anism, relies on the two-fluid effects that become important when the reconnection layer 



is exceeded 
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becomes so thin that ion and electron motions decouple from each other 
In reality, as the layer gets thinner, the Hall term becomes more and more important, but 
it is possible that the condition for anomalous resistivity is reached first. If this happens, 
then the enhanced collision rate due to fluctuations allows electrons to flow across the field 
lines instead of having to flow rapidly along them to preserve charge neutrality (see be- 
low). This will weaken or remove the Hall effect on reconnection. On the other hand, it 
may be that a stable two-fluid flow pattern is established first, before the instabilities that 
lead to anomalous resistivity are triggered. Then reconnection would proceed by the Hall 
mechanism. 

In the present study we focus on the Hall-magnetohydrodynamic (Hall-MHD) regime of 
reconnection. More specifically, our main objective is to understand the physical origin of 
the quadrupole pattern of the out-of-plane magnetic field that arises inside the reconnection 
region in this regime. We shall call it "the quadrupole field" for short. The quadrupole 
field is widely accepted as one of the most important signatures of the two-fluid effects 
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in the reconnection process. The reason for this is the following. Consider the simplest 
two-dimensional (2D) reconnection configuration (displayed in Fig. HJi), with z marking the 
ignorable direction, and x and y forming the so-called reconnection plane. Let us assume 
that there is no guide field, that is, assume that the magnetic field lines above and below the 
reconnection layer lie exactly in the reconnection plane ( "null-helicity" reconnection). Then, 
in simple resistive magnetohydrodynamics (MHD) there is no mechanism that would produce 
an out-of-plane (z) component of the magnetic field anywhere in the reconnection region; 
this is basically a consequence of the symmetries inherent to the resistive MHD equations. 
However, when electron and ion flows decouple from each other, that is when two-fluid effects 
become important, those symmetries are no longer present, since charge carriers of different 
sign now move differently. As a result, an out-of-plane component of the magnetic field may 
develop somewhere inside the layer. Thus, the emergence of this field is a tell-tale sign of 
the transition from resistive to two-fluid (e.g., Hall) regime of reconnection. For example, 
imagine a situation where one studies a reconnection process in a lab with limited diagnostic 
capabilities, say, with only magnetic probes but with no ability to measure plasma densities, 
temperatures, velocities, etc. Then one is not able to determine the ion skin depth and ion 
Larmor radius to compare them with the measured reconnection layer thickness. However, 
a mere detection of the ^-component of the magnetic field in the layer will immediately and 
unambiguously reveal that one deals with a two-fluid reconnection regime. 

The presence of an out-of-plane magnetic field with a quadrupole structure in the contex ; 



of the Hall-MHD regime of col 
(see also the work of Terasawa 



isionless reconnection was first suggested by Sonnerup 
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18j in the context of tearing instability in the Earth magne- 
totail). Since then, the quadru pole field has been observed in many numerical simulations 



of collisionless reconnection 
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22j. It has also been detected in space 



in-situ measurements by the Polar and Cluster spacecraft flying through the Earth magne- 



totail and the magnetopause 
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. The quadrupole magnetic field pattern has long 



evaded direct experimental detection in laboratory plasma experiments but just recently 



has finally been confirmed 



in the Magnetic Reconnection Experiment (MRX) 
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similar result has recently been reported in the Swarthmore Spheromak Experiment (SSX) 
in a somewhat different type of a neutral sheet (27]. These experimental detections raise the 
need for a better theoretical understanding of the generation mechanism for the quadrupole 
field, which in our view is still lacking despite the great wealth of numerical data. 

Our paper is structured as follows. In Sec. |Tl]we propose a basic physical explanation of 
how an out-of-plane magnetic field is generated in an electron MHD (eMHD) reconnection 
layer and why it inevitably has a quadrupole pattern. We start this section by discussing 
our physical assumptions (Sec. Ill A|) ; then we present a simple physical description of the 
mechanism by which the quadrupole field is produced fSec. IIIBj) ; and, finally, we illustrate 
our ideas by an analytical calculation of the toroidal (e.g., out-of-plane) field in the X-point 
configuration (Sec. Ill C|) . In Sec. |TTT| we step back from our specific example of Sec. Ill CI 
and derive some general results pertinent to stationary incompressible ideal eMHD in 2.5 
dimensions. Thus, in Sec. 1111 A1 we establish proportionality between three important quan- 
tities: the volume per poloidal flux computed along a field line, the poloidal (e.g., in-plane) 
electron stream function, and the electron contribution to the toroidal magnetic field. In 
particular, we show that, as long as ion currents are neglected in Ampere's law and the 
electrons are magnetized, the toroidal magnetic field is constant along electron streamlines; 
correspondingly, the toroidal electron velocity has to be constant along poloidal magnetic 
field lines. This means at least that inside the reconnection layer, at scales smaller than 
the ion inertial scale (but outside the inner electron dissipation region), one cannot invoke 



the usual explanation for the quadrupole field as being created by differential stretching of 
the poloidal field lines due to a non-uniform toroidal electron flow. It then follows that the 
toroidal field has to be generated in the transition region in the outskirts of the reconnection 
layer, where the ion-current contribution is still important. To study this process, we con- 
sider the three-dimensional shape of the field lines and show that the toroidal separation Az 
between a given fluid element on a line and the tip of the line at the x = plane is related 
to the volume per flux integral (Sec. IIIIB|) . This enables us to calculate the toroidal electron 
velocity (Sec. 1111 D|) and hence estimate when (i.e., how close to the separatrix) electron 
inertia becomes important in the generalized Ohm law (Sec. IIII E|) . Finally, we argue that 
the toroidal velocity of the field lines should be attributed to an E x B drift of electrons; 
this requires the presence of a bipolar poloidal electric field which we compute in Sec. IIII Fl 
This electric field is also an important signature of the two-fluid effects in reconnection; in 
particular, it is responsible for accelerating ions into the reconnection layer, resulting in an 
effective ion heating. We summarize our work in Sec. IIVI 

II. HOW IS THE QUADRUPOLE FIELD GENERATED IN THE RECONNEC- 
TION REGION? 

A. Physical Assumptions 

First, let us discuss the physical assumptions that we adopt in this paper. These assump- 
tions are aimed at making the problem tractable while still realistic and complex enough to 
provide a useful physical picture of a reconnecting current layer in the Hall-MHD regime. 
While doing this, we pay special attention to the conditions relevant to the MRX experiment. 

We will be mostly interested in the inner structure of the reconnection layer at scales (in 
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the direction across the layer) smaller than the ion collisionless skin depth defined as 



Provided that there is some ion heating available, and in the absence of a strong guiding 
field, the ion gyro-radius, p iy is comparable to dj in this region, and so ions can be regarded 
as unmagnetized. Thus, their motion is not strongly affected by the small-scale magnetic 
structures that characterize the inner part of the reconnection layer considered in this paper. 
On these small scales, the motion of ions is slow and smooth; the ion density then cannot 
develop structure on these scales. We shall therefore treat ions as providing a neutralizing 
background, which, for simplicity, we shall take to be uniform. Also, for the most part, 
we shall assume them to be motionless, that is we shall neglect the ion contribution to the 
electric current. However, as we will show, the poloidal ion current in the outer region of 
the layer actually plays an important role in the generation of the quadrupole field. 

In contrast to ions, electrons have very small gyro-radii and are well magnetized every- 
where except in a small vicinity of the X-point. Thus, it is appropriate to use the framework 
of Hall MHD (or electron MHD) in most of the region under consideration. This framework 
is characterized as a two-fluid approach where the magnetic field is frozen into electrons 
but not into ions. An equivalent formulation is to use a generalized Ohm law (i.e., electron 
equation of motion) that includes the Hall term. On the other hand, since we are interested 
in scales that are much larger than the size of the inner electron diffusion region, we shall, in 
our analysis of the generalized Ohm law, neglect both the electron inertia term (although we 
shall estimate its contribution in Sec. IIIIE|) and the resistive term that arises due to normal 
particle-particle collisions (i.e., classical Spitzer resistivity) or due to wave-particle collisions 
(anomalous resistivity). At the same time, we shall include the electron pressure gradient 
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term, assuming, however, that the electron pressure tensor is isotropic. This assumption is 
justified if the system is not entirely collisionless. That is, we assume that collisions are rare 
enough for collisional resistivity to be negligible, but, at the same time, frequent enough to 
restore the electron pressure isotropy throughout most of the reconnection region. This is 
in fact consistent with the conditions encountered in the MRX experiment, where collisions 
are always present at some level 26 1. 

Next, due to the charge neutrality condition (valid provided that the scales under consid- 
eration are still much larger than the Debye length), the electron density has to be equal to 
that of the ions. Since we assume that latter to be uniform, we require the electron density 
to be also uniform and hence the electron flow to be incompressible. 

Finally, we assume that the reconnection layer is in a quasi-steady state, that it has a 
translational symmetry in one (z) direction, and that there is no guide magnetic field (that 
is no externally imposed toroidal magnetic field). 

Thus, from the above considerations, the set of physical assumptions can be summarized 
as ideal incompressible 2.5-D steady-state electron MHD without a guide field 28 1. 

B. A Simple Physical Picture of the Quadrupole Field Generation 

We first describe the basic physical picture of how the quadrupole out-of-plane magnetic 
field naturally arises in electron MHD. Consider an incoming flux tube as it moves deeper and 
deeper into the (ion-scale) reconnection region toward the X-point (Fig. [T^). The poloidal 
magnetic field in the central part of the tube near x = has to decrease, and hence the 
volume of this central part has to expand. Since electrons are tightly coupled to magnetic 
field lines, this expansion would lead to a drop in electron density. However, the ions are 
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not magnetized and their density does not decrease. Therefore, since almost perfect charge 
neutrality is to be maintained, a very small poloidal electric field arises and it immediately 
pulls the electrons along the field lines inward from the outer parts of the flux tube into 
this central region. Owing to the very large mobility of electrons along the field (inversely 
proportional to m e ) , this parallel electric field is negligibly small. 

As a result, we get a strong inflow of electrons along the poloidal magnetic field in the 
upstream region (Fig. [TJ)). This inflow rapidly accelerates as the field line approaches the 
separatrix, because of the rapidly increasing rate of flux-tube expansion near the X-point. 
In the downstream region, the direction of the electron flow reverses: as a newly reconnected 
field line moves away from the X-point, the volume of its central part decreases and so the 
electrons are squeezed out and flow rapidly outward along the field (Fig. Ct). As the field 
line moves further away, this outflow gradually decelerates. The resulting overall picture of 
the electron flow is shown in Fig. once again, the main feature is the rapid inflow of 
electrons just above the separatrix followed by a rapid outflow just below the separatrix. 

This pattern of electron motion plays an important role in eMHD reconnection, since 
there is a poloidal electric current associated with the flow of electrons. By Ampere's law, 
this current generates a quadrupole toroidal magnetic field concentrated along the separatrix 
(see Fig. |TJi). This is our picture for the origin of the quadrupole field. 

The orientation of this field is always such that the toroidal field in the upper right and 
lower left quadrants is directed away from the viewer, whereas the toroidal field in the lower 
right and the upper left quadrants is directed towards the viewer. It is interesting to note 
that this orientation is universal, i.e., independent of the direction of the poloidal field. 
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C. An Analytical Example: a Simple X-point Configuration 

To illustrate this mechanism, we present a very simple calculation of the toroidal field 
based on the simplest possible poloidal field configuration relevant to the reconnection prob- 
lem. This configuration is of course the X-point configuration that we now describe. 

Consider the central part of a reconnecting current layer. Let L be the half-width and 
(5 < L be the half-thickness of the layer. Let us choose a Cartesian coordinate system 
(x, y, z) with x being the direction along the layer, y across the layer, and z in the ignorable 
direction. We shall refer to the plane of reconnecting field (i.e., the xy plane) as the poloidal 
plane and the z direction as the toroidal direction. We set the origin x = = y exactly at 
the X-point and assume mirror symmetry with respect to the xz and yz planes as well as 
the translational symmetry in the z direction (see Fig. |2J). We also assume steady state. 

We are interested in a small vicinity of the X-point, which means that we consider loca- 
tions with i < i and y <C 5. In this region the poloidal magnetic field can be generically 
represented by a simple X-point configuration. In terms of the poloidal flux function \l/(x, y, t) 
this can be written in the appropriate gauge as 

»(,,**) = -^^(g-g), (2, 

with B poi = W x z. 

This expression serves as the definition of the scales L and 5. The first term in this 
expression is just the instantaneous value of the flux at the origin (which we can define 
as the flux that has reconnected since t = 0). As reconnection proceeds, it increases at a 
constant rate equal to —cE z = \cE z \ > 0, where E z < is the toroidal electric field (which is 
uniform in steady state reconnection). The quantity Bq represents the reconnecting magnetic 
field just outside the layer: B x (x — 0,y — 6) — B . The poloidal magnetic field components 
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corresponding to this flux function are: 



B x = d y m = B | , (3) 
o 



B y = -d x ^f = B ~. (4) 



For convenience, we introduce dimensionless variables by rescaling x, y, B, and E z as 

x y 

B = * - cE z 

Then equation (J2J) can be written as 

—2 —2 

y> = y + Et = *(x, y) - f (0, 0) = V - - X — . (6) 
Correspondingly, the shape of a given field line ^ is given (in the two upper quadrants) 

by 



y(x,V) = v / 2^ / + x 2 . (7) 

Now we want to calculate the motion of the electron fluid. It is completely determined 
by two conditions: flux-freezing in the poloidal plane (which is not spoiled by the pressure 
gradient term in the generalized Ohm law, as we shall discuss later) and incompressibility. In 
order to get an explicit expression for the electron velocity, let us consider the trajectory X(t), 
Y(t) of an electron fluid element. As it moves through the layer, the given fluid element 
always stays on a field line with constant thus, as one follows its motion, = ^ + Et 
varies with time. Correspondingly, the trajectory of the element has to satisfy 



Y(t) = y/2(* + Et) +X 2 {t). (8) 

Next, the incompressibility condition implies that the volume per unit flux following an 
electron fluid element, V(X, = V[X(t), ^ + Et], has to be conserved. Here, V(X, is 

10 



measured along the line from the y-axis (x = 0) in the case of a field line in the upstream 
region (and from the x-axis, y = 0, in the case of a field line in the downstream region), up 
to the fluid element under consideration. For example, in the upstream region we thus have 



l(X) 



V[X(t),V] 



dl 



o 



B 



pol I 



dx 

\t=const J \^x\ 



= const . (9) 

>&=const 



Using the expression (J3J) for B x and (J7J) for the field line shape y(x, we get 
x(t) 



v\x{t)M = T J 

o 



L f dx 

o 



= log 

*=const -Do 



X(t) + JX 2 (t) + 2W(t) 



const . (10) 



Hence, the trajectory is given by (X(t) + Jx 2 {t) + 2^')/v / 2$ 7 = const, that i 



is 



X(t)/ J2^/'(t) = const. Using (jHJ), we also get a similar expression for Y{t). Thus, 



X(X ,t) = £(*o,*W2|* + £*l, (11) 



Y(X ,t) = V (X ,*)y/2\* + Et\, (12) 

where 77 = V£ 2 ± 1, and where we take "+" in the upstream region and "-" in the downstream 
region. The constant parameters £ and r] represent the initial position of the electron fluid 
element at t — 0. 

Differentiating these expressions with respect to time and using d^'/dt = E, we obtain 
the electron velocity field at any point (x, y): 

-(e)/- , x E _ \E\ 

Vs \x,y) = - -._ , = -x— — , (13) 



'Ti-r.y) = f — f— = -y -13- . (14) 



2 ty'(x, y) y 2 — x 2 

y E \E\ 
2 ^/'(x, y) ^ y 2 — x 2 

Notice that vf^ jv^ = y/x, so electrons flow along purely radial lines, inward in the 
upstream region and then outward in the downstream region. Even though all the streamlines 
converge to the X-point in the upstream region (and fan out of the X-point in the downstream 
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region) the motion is incompressible since the magnitude of the velocity diverges near the 
origin. Also note that if the above velocity field holds, electrons never actually cross the 
separatrix; they all go through the X-point. This is of course an artifact of our ideal eMHD 
assumption. 

This poloidal electron velocity field results in a poloidal electric current: 

J x e) = ~ en e v e x = - en e Lv% , (15) 
jy e) = - en e Vy = - en e 5v e y . (16) 

This current in turn produces a toroidal magnetic field. According to Ampere's law we have 

d S M = — W = ^ L\E\ I (17) 

y z c 3x c 1 l 2V(x,y)> K n 



c Jy c 1 1 2V'{x,y) ' v ' 

We can thus compute the resulting toroidal field by integrating either of these equations 
while taking into account that B z (0, y) = = B z (x, 0) because of symmetry. In the upstream 
region it is convenient to integrate d x Bf'(x,y) with respect to x at constant y starting with 
x = 0. (In the downstream region, it is convenient to do the opposite.) Thus we can write: 



X X 



B z (x, y>x) = J Ld x B z dx = -- ^ SL \E\y J . (19) 

o o ^ ,y) 

Using the expression (JBJ) for ty'(x,y), we get 



fl,(x,jO = -is Qlog|^?|, (20) 
2 y — x 

where we define a dimensionless constant 

Q= ^6L\E\ = ^\E\. (21) 
cB di V A 

Expression (|20|) is actually valid in both the upstream and downstream regions. We note 
that a very similar expression was obtained, in the eMHD framework, in Ref. 
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We can rewrite the coefficient Q in a different form by expressing the reconnection electric 
field E in terms of the other parameters of the reconnection region. Indeed, let v TCC be 
the reconnection velocity, v Tec = —v y (x = 0,y ^> 5). Then \E Z \ = v rec B /c and hence 
\E\ = v rec /5. Then, we get 

= ^ = 4" (22) 

di V A a,, V A 

where u is the velocity of the flow out of the layer, and where we define a new dimensionless 
parameter C: 

C=^ = 0(1). (23) 
ou 

Because of the condition of overall mass conservation we expect C to be of order unity. 

Usually, one expects the thickness 5 of a Hall-MHD reconnection region to be comparable 
to the ion collisionless skin depth di and the outflow velocity u to be of order Va- Thus, 
equation tells us that the proportionality coefficient Q is, generally speaking, expected 
to be of order unity. In practice, however, Q may significantly deviate from unity for any 
specific physical system. For example, in the MRX experiment one often encounters 5 ~ di/3 
and u < Va, so Q can be smaller than 1/3. 

The reader should be warned that formula (J2(J|) only applies to our specific example 
for the poloidal field, eq. (j2J), and may not be applicable to various configurations realized 
in some numerical simulations and in the MRX. We chose this specific example because 
of its simplicity and clarity and we leave more complicated poloidal field structures for a 
future study. We believe that the most physically-relevant among these other structures is 
a configuration with an inner electron current sheet (electron dissipation region) of a finite 
width in the x-direction. 

As we see from equation (J2"U|) . in the simple X-point configuration considered in this 
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section, the assumptions of ideal eMHD lead to a logarithmic divergence of B z at the sepa- 
ratrix y = x. Later, in Sec. lIIIEl we shall discuss how this singularity is removed by including 
finite electron inertia. We also see that B z in our solution is constant along straight radial 
rays y = const • x. In reality we expect electron inertia and other non-ideal effects to in- 
tervene and break this idealized picture near the separatrix. However, we believe that the 
main tendency for the constant-.B z contours to be strongly elongated along the separatrix 



will survive. In fact, t 
numerical simulations 

data 



lis overall behavior is in a very good agreement with the results of 



15 



16 



19 



20, 
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221 ] . and is also consistent with the experimental 



III. STATIONARY IDEAL INCOMPRESSIBLE EMHD IN 2.5 DIMENSIONS: 
GENERAL RESULTS 



In this section we step back from the particular example of the previous section and 
derive several general results that are valid in steady-state 2D incompressible eMHD for an 
arbitrary poloidal field structure. 



A. General relationships between toroidal magnetic field, electron stream function 
and the volume per flux in eMHD 

First we introduce three important functions: the volume-per-flux integral V(x, the 
electron stream function <3> e and the (electron contribution to) the toroidal magnetic field B z . 
We derive important relationships between these functions and discuss their implications. 

The volume-per-flux integral V(x, is defined (in the upstream region), as in the last 
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section, by 



dlp D \ 

I -Spoil 



(24) 

1 1 

where, in the last expression, the integration is performed along a given poloidal line, 
from the y-axis x = to the given point (x, A similar expression can be defined in the 
downstream region. 

The electron stream function $ e (x, y) is defined, for an incompressible flow, V • v* 6 ) = 0, 
by the poloidal electron velocity as 

vS = [V x (<M)] = [V$ e x z) . (25) 

For definiteness, we choose the streamline corresponding to zero $ e to coincide with the line 
from which we count the volume-per-flux, i.e., with the ?/-axis: $ e (0, y) = 0. 

With these definitions we now show that in steady state ideal incompressible electron 
MHD the two functions are just proportional to each other: 

$ e (x,y) = -cE x V(x,y). (26) 

where E z is the uniform electric field in the ^-direction. 

To prove this, let us consider the variation of $ e along the poloidal magnetic field and show 
that it is proportional to that of V. Consider two points lying close to each other on the same 
poloidal field line and let Al po \ be the infinitesimal separation between these two points 
along the poloidal field. From equation (|24p. the difference between the volume-per-flux of 
these two points is 

AV = ^ . (27) 
Bpoi 

On the other hand, the difference between the values of the electron stream function at 
these two points can be expressed in terms of the perpendicular component of the poloidal 
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electron velocity 

v£L = [V$ e xf] i = ^ [b pol X |] = - v * . ( 28 ) 



where b po i is the unit vector along the poloidal magnetic field. On the other hand, ± 
can be deduced from the generalized Ohm law. Neglecting resistive and inertial terms and 
taking into account that d z p e = 0, the toroidal component of this law can be written as 

cE z z+[v^ x xB pol } = 0. (29) 

By taking the vector product of this equation with B po i, we get 
( e ) E z z x B pol E? r „ , E 



v P ol,± - c - c ~ "poij - ~B2 

Comparing equations (|2*8j) and (|5U|) we immediately see that 

A$ e = - — - AZpoi , (31) 

-Dpol 

which, together with equation (|27|). gives 

A$ e = - cEgAF . (32) 

Because we choose to count the volume-per-flux from the $ e = electron streamline, 
then, summing this equation along each field line, we get equation ([26)1. as desired. But we 
can actually proceed more directly. Let us consider the variations of $ e and V across the 
poloidal field. For this, consider a single poloidal field line \I/ at two neighboring moments 
of time, t\ and t 2 = t% + At. Consider a certain point A\ on the field line at t — t\ and see 
it E x B-drift with the field line to a new location A2 at t = t 2 - The corresponding change 
in the volume-per-flux is due to the plasma that has flowed in along the poloidal field past 
this point during the time At: 

AV = V(A 2 ) - V{A X ) = - . (33) 

-Dpol 
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[z x b pol ] = c -f - Vtf . (30) 



At the same time, this parallel inflow of the plasma results in a change in $ e between 
points A\ and A 2 : 

A$ e = $ e (A 2 ) - $ e (A 1 ) = v poU (As poU • [z x bp i]) , (34) 

where the poloidal displacement vector As poliJ _ in the direction perpendicular to the poloidal 
magnetic field is given by the E x B drift: 

cE z 

As po1i _l = v_L iPol At = At — — [z x bpoj . (35) 



Thus, we get 



Bpoi 



cE 

A$ e = u pol) || At — ± , (36) 

-Dpol 



and comparing this result with equation we again see that 

A$ e = — cE z AV . (37) 

A similar derivation also holds in the downstream region. 

Thus we have shown that the variation of the electron stream function $ e in both parallel 
and perpendicular directions is equal to —cE z times the corresponding variation in the 
volume per flux integral V. Using the convention of counting both $ e and V starting from 
the y-axis, we again arrive at the relationship (|26j) . 

The second important relationship is the proportionality between the electron stream 
function $ e and the electron contribution to the toroidal field B z . This well-known re- 
lationship follows immediately from Ampere's law and the reflection symmetry conditions 
[B z (x = 0,y) = = $ e (x = 0, y) upstream and B z (x,y = 0) = = & e (x, y = 0) down- 
stream]. It reads: 

B z (x,y) = -D<S> e (x,y), (38) 
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with the coefficient D given by 



„ 4iin e e \/4:TTp Bn , , 

D = e - = ^— H - = — °- . 39 

c di diV A 

(Here Bq is an arbitrary normalization field used in the definition of Va-) 

Note that the coefficient D defined by equation is constant for the case of uniform 

density considered here. If the density were not uniform, we would get a similar result B z ~ 

$ e , if the electron density is incorporated into $ e , i.e., if <3> e is defined by UevQ = [V x (3> e z)]. 

Similarly, equation would be valid if V is understood as the number of electrons per 

unit flux, instead of the volume-per-flux. 

Combining this equations (|38|l and we immediately see that 

B z = cDE z V. (40) 

This result is important because it shows how, in eMHD, one can immediately determine B z 
once the poloidal field structure is known, without having to solve any partial differential 
equations! One just has to compute the volume-per-flux integral given by equation 
upstream and the corresponding expression downstream, and the toroidal field will be just 
the constant cDE z times it. 

Already by itself, the simple relationship (J38)) is important, because it means that, ne- 
glecting ion currents, the toroidal field is constant along the poloidal electron streamlines; 
its value is simply transported in space by the poloidal electron flow: 

(vg ■ V)B Z = - DB [V$ e xi].V$ e E0, (41) 

where is the total poloidal electron velocity including the parallel flow. 

This result suggests that the toroidal magnetic field cannot be created locally, in the 
inner part of the reconnection region (where ion current is unimportant); instead, it has to 



be brought into this region by the convergent electron flow. The toroidal field thus has to 
be generated in the outer parts of the layer, where the ion-current contribution to B z is not 
negligible. We discuss this generation process in Sec. 1111 CI 

In addition, as long as electrons are completely frozen into the magnetic field, the evo- 
lution equation for the toroidal magnetic field, i.e., the toroidal component of the magnetic 
induction equation, tells us that 



v. 



(<0 
pol 



• VB Z = B pol ■ Vvi e) , (42) 



in a steady state. Therefore, since the left-hand side (LHS) of this equality is zero, as we have 
just shown, the right-hand side (RHS) is also zero. That is, the toroidal electron velocity and 
hence the toroidal current density j z = — en e v^ are uniform along poloidal field lines as long 
as one can neglect the ion current. Since from Ampere's law the toroidal current density is 
simply proportional to the Laplacian of we see that in order to get a consistent solution, 
one cannot pick the poloidal flux function arbitrarily; one has to impose the condition that 
V 2 \I/ = F(ty). For example, the simple X-point poloidal field configuration considered in 
Sec. Ill CI trivially satisfies this condition with Ffif) = const. 

The above observation also means that within the pure eMHD framework with no ion 
currents one cannot really apply the well-known conventional explanation, first introduced 



in Ref . 14j , of how the quadrupole field is generated. By itself, this argument does not rely 
on neglecting ion currents. Instead it relies on the fact that magnetic field is completely 
frozen into the electron fluid. The toroidal field is then viewed as being produced from the 
poloidal magnetic field as a result of the differential stretching in the toroidal direction by 
the electron flow. This argument thus approaches the generation of the quadrupole field 
from a different angle: it presents the point of view of the ideal eMHD Ohm's law (with 
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the Hall term), instead of using Ampere's law. It basically goes like this: as a field line is 
advected into the reconnection layer, the electrons on it start to move toroidally (to carry 
some of the reconnection current) and they do it differentially, moving faster on the central 
piece of the field line. Since the magnetic field is frozen into the electron fluid, the field line 
bends out of the reconnection plane, resulting in the quadrupole pattern of the toroidal field. 
This line of thought is actually quantified by equation (|%T|) . That is, the toroidal field is 
created by the stretching due to the non-uniformity of along a poloidal field line (RHS) 
and is advected with the poloidal electron flow (LHS). But, as we have just seen, as long as 
one neglects ion currents, this equation becomes simply = 0. Thus, in order to understand 
toroidal field generation, one first needs to take ion currents into account (see Sec. IIII CJI for 
more discussion). 

B. The Shape of Field Lines in the xz Plane 

A complimentary way to look at the problem of the toroidal field generation is to analyze 
the shape of a field line projected on the xz plane, and see how it changes as the field lines 
move deeper into the layer. Let this shape be represented by the function z(x, which is 
given by 

dz . B z , . 

1* = -^- (43) 



dx B 



X 



By integrating this along a field line we obtain: 



/dx 
B Z (X,V) , Tf J *=const- (44) 



Using equations (jlUj) and ({21$ , we then have 



V 2 (x \f0 

Az{x, *) = cDE z { ' ' . (45) 
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Note that the volume per flux is conserved by the motion of the electron fluid; this means 
that the x-position of an electron fluid element that stays on some constant- 1 !/ field line 
changes with time in such a way as to keep V(x, constant. Thus, if one follows a specific 
fluid element on a given moving field line one finds that the toroidal distance Az(x, 
between this element and the point where the field line intersects the x = plane does not 
change with time. On the other hand, as we showed in Sec. Ill Bl the fluid element moves along 
the poloidal magnetic field towards the y axis, and so its x-coordinate decreases. Therefore, 
the shape of the field line, which can be characterized by the function x(Az), changes with 
time. Interestingly, this change is not due to the differential toroidal stretching, as it is 
usually assumed, but is simply due to the fact that the two mirror-symmetric parts of the 
line are squeezed together by the converging poloidal flow. 

C. The Role of Ion Currents in the Generation of the Quadrupole Field 

As we saw in Sec. IIII A| when the current due to the ions is neglected, the toroidal field 
is conserved along the electron streamlines and hence cannot be locally generated in the 
inner part of the reconnection layer. This indicates that, in order to explain how and where 
the toroidal field is generated, one has to bring the ions back into the picture. Deep inside 
the reconnection region, at x <C L and y <C 5 ~ dj, the poloidal ion currents are indeed 
negligible and the above picture applies. On the other hand, in the upstream region well 
outside of the reconnection layer (i.e., for y ^> 5), ideal one-fluid MHD works well. In 
this region the electrons do move poloidally towards the reconnection layer with the E x B 
velocity and the associated electron current would generate the toroidal field; however, the 
ions also move happily in the same direction and with the same speed. As a result, the 
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ion-current contribution to the toroidal field exactly cancels that of the electrons. Thus, the 
net toroidal field is zero in this region. From this we see that, in order to understand where 
the quadrupole toroidal field comes from, one has to look at the outskirts of the reconnection 
layer, where the ions become partly decoupled from the electrons, so that < < 

We can discuss the toroidal field generation from a different point of view, in terms of the 
shape z(x, of a given field line \I> as it is carried into the current layer by the electron flow. 
Far upstream, this field line lies entirely in the reconnection (x, y) plane, but as it moves into 
the reconnection region, it gradually starts to bend out of this plane. The toroidal electron 
velocity can be non-uniform along the line only in this transition region of non-zero ion 
current. Correspondingly, toroidal field is produced inside this region; subsequently, deeper 
inside the reconnection layer, the toroidal electron velocity becomes uniform along the line 
and hence the toroidal elongation freezes. Any further lengthening of the field line in the 
toroidal direction can be directly attributed to the "injection" of new segments of the field 
line in the transition region. 

D. Toroidal Electron Velocity 

To illustrate this picture, let us consider an extremely simplified model where the transi- 
tion region is a razor-thin line y — 5. In this example, the electrons and ions move together 
above y = 5 and so B z (x, y > S) = 0. Below this sharp boundary, we shall regard the ions as 
poloidally motionless, so that j®i(y < 5) = and hence the pure eMHD picture developed 
in the preceding sections applies. In addition, in this and in the next section we shall, for 
simplicity, neglect the toroidal component of the diamagnetic electron flow that results from 
poloidal electron pressure gradient. 
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Now let us consider a given field line as we follow its motion through the layer, the 
magnetic flux ty' between the separatrix and the given field line changes linearly in time 
according to = \P + cE z t. Denote the x-coordinate of the point where this field line 

intersects the boundary y = 5 by ^[^'(i)] = x s{^ + cE z t). For example, in the simple 
X-point configuration considered in Sec. Ill CI we have Xg — [1 — 2^f'(t)] 1 ^ 2 . Next, because 
there is no toroidal field above y = 5, we can set z = everywhere along this boundary, i.e., 
z[xg(^'), = 0. Then, using equation (f4~HJ) . we can express the toroidal coordinate of any 
fluid element {X(t), f) on a given line \1/ as 



i[X(t), *] = Az[X(t), *] - Az[x 5 (*'), *] = cDE, 



^ 2 [X(t),^]-^ 2 (vI/') 



(46) 



where 



V s (*')=V[x s (*'(t)),*]. 



(47) 



In this section we are interested in the toroidal electron velocity, so let us see how z(X, \I/) 
changes with time following an electron fluid element. To do this, differentiate equation (|4fjjl 
with respect to time. When doing this we have to take into account that in ideal incompress- 
ible electron MHD the motion X(t) of a given fluid element is constrained by the condition 
that V[X(i),\l/] remains constant. Then we have 



vim 



c 2 DE 2 



d 



(4£ 



since d^' /dt = cE z . Thus, the velocity is proportional to the flux derivative of the square of 
a flux tube's entire volume up to the boundary y = 5. 

One sees that the toroidal velocity is constant along field lines but, in general, varies from 
line to line. In particular, the volume-per-flux Vg(^f') grows rapidly near the separatrix and so 
v^f 1 becomes very large there. This appears to be inconsistent, for instance, with the simple 
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X-point configuration considered in Sec. Ill C| indeed, the particular form of the poloidal 
flux function in that example corresponded to a flat toroidal current profile, j z = const, and 
hence v[ e ^ = const. The way to resolve this discrepancy is to note that the sharp rise in the 
toroidal current density that corresponds to equation (|48|). leads to only a relatively small 
change in the poloidal field structure. Moreover, this change is actually consistent with that 
expected from the back-reaction of the toroidal magnetic field pressure. This back-reaction 
arises because, as one approaches the separatrix, the toroidal field increases sharply and 
starts to play an important dynamical role. In particular, it modifies the poloidal field 
structure through the vertical pressure balance condition; the poloidal field decreases near 
the separatrix and this leads (by Ampere's law) to an additional electric current, strongly 
concentrated near the separatrix. We can estimate this additional electric current as follows. 
Let us write the vertical pressure balance as 



where we have assumed that the total plasma pressure P has a parabolic profile: 8tiP = 
Bq(1 — y 2 ). If we neglect the toroidal field pressure term in this equation, we then recover 
our original poloidal field profile B po \ « B x = Boy, which corresponds to uniform toroidal 
current. Note, however, that even if B x itself is small, its rate of change may become 
important near the separatrix, so that the corresponding small but rapid change in B x 
results in a large contribution to the toroidal current. Indeed, differentiating equation (J4U| 
with respect to y, we get 



BL + Bl = Bl - 8nP = Blf , 



(49) 



dBl 



pol 



dB x 



BJj z = 2B 2 y 



dBl 



(50) 



dy 



2B X 



dy 



dy 
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Then, using equation @ and (pfUj) . we find 

'V 2 (x,^') 



3z ~ 4tt 



T~ cDE *d¥ 



n P ec 2 DE 2 



d 



V 2 {x^') 



(51) 



and so, assuming that the additional toroidal current is predominantly carried by electrons, 



vP~-(?DEl 



d 



v 2 (x, m 



(52) 



- an expression that is very similar, although not quite the same, as equation 
Finally, for reference, let us give expressions for Vs(^f') and that correspond to the 
simple X-point configuration considered in Sec. Ill CI First, according to equation (JTUJ), we 
have 



L 



Br 



log 



1 + Vl - 2^' 



= log 


1 


2B 





so that (d/dW) V S (V') = -(L/2B ty'x s (y') 
Then, from equation (|48|) we obtain 



= - B 5D\E\ 2 V s (*') V s (*' 



d 



L 2 5D \E\ 2 



dV ' 2B 2V'x s (V') 
Finally, using definition (|21jl . we can write this as 



l0£ 



i-x s m 



LQ\E\ 



7T log 



i-^(*0 



(53) 



(54) 



(55) 



E. Finite Electron-Inertia Effects 

We can use the above formula for v^f 1 to estimate when the electron inertial term stops 
being negligible in the toroidal component of the electron equation of motion. From that 
moment on, the finite electron inertia will be large enough to balance part of the toroidal 
electric field, and thus, the electrons will no longer have a pure E x B velocity and will no 



25 



longer follow the field lines exactly. The inertial term (for a single electron) can be written 

as 

m e (v • V) « m e (v pol , ± • V) v& = m e (v poU • W) , (56) 

where we take into account that v}?' is constant along field lines and so is a function of 
only. Using expression (|50|l for v po i,±, we get 

m e (v • V) vf « m e c£ z . (57) 

Then, using (pfSj). we get 



m e (v • V) « - m e c 3 E 3 z D 



d 2 rv 5 2 {V) 



(58) 



d(#') : 

We can no longer neglect electron inertia when this term becomes comparable to the 
toroidal electric force on an electron, — eE z . We estimate this to happen for values of ^f' of 
order of which is obtained as the solution of the equation 

= J ^4w-2 = v^— 4* ■ ( 59 ) 

m e cD c^E* 27rn e m e crE* 

We can apply this estimate to our simple X-point example, for which Vg(^f') is given 
by (|53p. Since we expect the electron inertia to become important only near the separatrix, 
\&' <C 1, we can approximately write Vg(^f') ~ —(L/2B Q ) log^', and then 

v z e) « - - QL\E\ ^f- . (60) 
2 2\l/ v ; 

The inertial term in the toroidal component of Ohm's law is then estimated, with the help 
of equation ()57|). as 

/f,,(e) i np . ^y/ 

m e (v • V) «W « m e £ « - m e QL|£| 2 ^ , (61) 
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where we recall that E < in our solution. After some manipulation we can write the 
condition on ty' as 



2% _ g de _ \rnl Lv iec _ ^ 
logvj/'J 5 V m * 5Va 

where we use equation ({22} and define d e = c/u pe . Thus, 



U,; 



(62) 



^~J-^log-i. 63 

V "It "^e 

When traced to the y axis, this critical field line corresponds to a distance 



~ — log — 5 64 

\rrii m e J 

from the X-point. (This is of order 5/4 for hydrogen plasma). It is essentially (apart 
from the logarithmic factor) of the same order as the distance at which electrons become 
demagnetized, i.e., comparable to the size of electron figure-eight and betatron orbits. 

F. The bipolar poloidal electric field 

Why do field lines move in the toroidal direction as they enter the layer? To answer this 
question, we need to consider the toroidal projection of the perpendicular (to the total mag- 
netic field) electron velocity, v± jZ . Let us locally introduce a rotated orthonormal coordinate 
system (x', y', z) where x' is the direction along the poloidal magnetic field and y' is the 
direction in the poloidal plane which is perpendicular to the poloidal magnetic field. Taking 
into account the electron pressure (which we assume isotropic) but neglecting the electron 
inertia, we can express v± z as 

(e ) [Eppi x B P oi] z [Vpoi(p e /n e e) x B po i] z 



B 2 B 2 
EyiBp i d y '(p e /n e e) Bpo\ 

- c ^h~ c - — & — • (65) 
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This is a sum of two drifts: the E x B-drift due to the poloidal electric field E y > and the 
diamagnetic drift due to the electron pressure gradient. In principle, as long as the electron 
pressure is isotropic, these two terms can be combined by noticing that in a steady state 
the poloidal electric field is electrostatic, E po i = — V02(^, y), and defining 2 = 02 — p e / n e e - 
Then, 

vf z = c ^faB** (66) 



However, an important point is that the diamagnetic drift is actually irrelevant, as far 
as the motion of field lines is concerned. In the presence of the pressure gradient, the field 
line velocity in fact differs from the electron perpendicular velocity and is given by just the 
E x B velocity. Its z-component is 

v B ,! = -c^. (67) 

Thus, the field lines move toroidally because of E y / that has a bipolar structure (see 
Fig. EJ). The above argument suggests that, instead of saying that electrons pull the field 
lines in the toroidal direction in two-fluid reconnection, it is, in a sense, better to say that it 
is the magnetic field lines that start moving toroidally and pull the electrons with them. The 
poloidal electric field can therefore be viewed, similarly to the quadrupole toroidal magnetic 
field, as an important signature of Hall reconnection. It has in fact been detected with with 
the Polar spacecraft in the magnetopause 



regions in the Earth magnetotail 
seen in numerical simulations 
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with the Cluster spacecraft in reconnection 

n 



and in the SSX experiment |27[; it has also been 
2l|. It is this electric field that pulls ions into the 
reconnection layer; as they move across the layer, they pick up the elecrostatic potential 
difference of the order of 5E y > . This potential difference is large enough to accelerate ions (in 
the ^-direction) up to about Alfven speed. As a result, the ion ^-distribution at the center 

28 



of the reconnection layer is well represented by two counter-streaming beams, which a gree s 



both with numerical particle simulations 



17 



201 ] and with spacecraft measurements 



24 1 . 



Effectively, this process can be interpreted as a strong ion heating, providing the pressure 
support for the layer. In addition, ion collisions (with particles or with waves) may quickly 
isotropize the ion distribution function, leading to a true ion heating. The quadrupole 
toroidal magnetic field also plays an important role in the poloidal ion motion; in particular, 
as the ions are accelerated into the layer, the Lorentz force due to this field bends ion 
trajectories in the ^-direction and thus leads to the ejection of ions out of the reconnection 
region. 

For reference, we give an expression for the bipolar electric field for our simple X-point 
configuration example. To derive this expression, we make use of the toroidal electron 
velocity vp, computed in Sec. 1111 151 

First, from the ^/-component of the ideal eMHD Ohm's law, we can write (neglecting 
electron pressure) 



cE y , = -vfB x , + v£>B. 



,(«0 



5poi + v poUl B z = -B pol {vf - \B Z ) 



(66 



where we express v x > = f po i,|| as \B po \ as it is done in Appendix |XJ In the case of the simple 
X-point configuration of Sec. Ill C[ we have at our disposal explicit expressions for all the 
ingredients that enter the above equation. Thus, B po \ is approximately equal to B x given 
by equation Q; v^f* is given by equation (j3SJ), A by (jA10|) . and B z by (j2Uj) . Putting it all 
together, we obtain 



y 



E, u i(x, y) = — — \E Z \ — —lr , 
y v ,yj 2 1 1 5 2^' \xs(^ 



77 log 



x log 



y + x 



y-x 



(69) 



We see that, because of the L/5 factor, this poloidal field can be considerably larger than 
the toroidal electric field. 
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IV. CONCLUSIONS 



In this paper we have investigated the structure of a reconnection layer in the Hall- 
MHD regime, in which electrons are well-magnetized inside the layer, whereas ions are 
not. Specifically, we have addressed the issue of how the quadrupole pattern of an out- 
of-plane (toroidal) magnetic field is generated inside a Hall-MHD reconnection region. This 
quadrupole pattern is commonly seen as an important feature of two-fluid physics that is 
at work in the reconnection process whenever the resistivity is small. It has been routinely 
observed both in numerical simulations and in space, and has recently been confirmed in 
a dedicated laboratory experiment |26j]. In our view, this quadrupole pattern arises most 
naturally via the following mechanism. 

Let us follow a flux tube as it enters the reconnection layer from the upstream region. As 
it moves deeper into the layer, the in-the-plane (poloidal) field in the central part of the tube 
weakens and so its cross-sectional area expands. This does not affect the ions very much. Let 
us assume that their density is constant throughout the inner part of the reconnection layer. 
Then, owing to charge neutrality, the electron density also has to be constant. Therefore, 
since the central part of the flux tube is expanding, electrons have to flow in into the layer 
along the poloidal magnetic field. Similarly, in the region downstream of the X-point, the 
flux tube is leaving the layer and so its cross-sectional area contracts. The electrons then are 
forced to flow along the poloidal field out of the layer. We thus obtain a circulating pattern 
of the electron current. In turn, it gives rise, through Ampere's law, to a toroidal magnetic 
field that automatically has a quadrupole structure. A more detailed qualitative description 
of this process is presented in Sec. Ill Bl 

We find that the most elegant and effective way to quantitatively analyze the behavior 
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of the system is in terms of the volume-per-fTux integral V(x, which has a nice property 
that it is determined entirely by the poloidal magnetic field structure, ^(x, y). We show that 
both the electron stream function <3> e and the toroidal magnetic field B z are just proportional 
to V . Thus, once the poloidal field structure is specified, the poloidal electron velocity and all 
three components of the magnetic field are easily determined just by computing one integral, 
i.e., without solving any partial differential equations. In particular, we find that, as long 
as poloidal ion currents are neglected, the toroidal magnetic field is constant along electron 
streamlines. This means that, within a pure eMHD framework, the toroidal field cannot be 
produced! Instead, it has to come from the outer regions of the reconnection layer, where 
ion currents are not negligible. 

We also find that the toroidal magnetic field is highly concentrated near the magnetic 
separatrix. We obtain explicit expressions for a simple X-point configuration and show that 
B z has a logarithmic singularity at the separatrix. In reality, of course, the vertical pressure 
balance condition would prevent the toroidal field from being larger that the outside poloidal 
magnetic field B . This means that the pressure associated with the toroidal magnetic field 
becomes dynamically important near the separatrix and hence the poloidal field structure 
must be such as to keep the toroidal field finite. In addition, the singularity at the separatrix 
is removed by electron inertia. In order to estimate the electron inertial term in the toroidal 
component of the generalized Ohm law, however, one needs to know the toroidal electron 
velocity. To determine it, we consider how the full three-dimensional shape of a field line 
changes with time as the field line is advected into the layer. From this we deduce the 
toroidal electron velocity and hence estimate how rapidly the electron inertial term grows 
near the separatrix. This enables us to estimate size of the region around the separatrix 
where the electron inertia is not negligible. 
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It should be remarked that, in spite of the fact that our calculation diverges near the 
separatrix, it is perfectly valid in the upstream and downstream regions away from the 
separatrix. In fact, all our integrations are carried out from the x and y axes towards the 
separatrix and do not cross it. 

Finally, we consider the well-known explanation of how the quadrupole toroidal field is 



produced (see, e.g., Ref. JJ|). This explanation invokes the differential stretching of poloidal 
field lines by a non-uniform electron flow in the toroidal direction. So a natural question to 
ask is: what makes the electrons move in the toroidal direction inside the layer? We argue 
that the toroidal electron velocity is in fact the sum of the E x B drift, associated with 
the bipolar poloidal electric field that points into the layer, and the diamagnetic drift due 
to electron pressure gradient. However, the latter does not lead to any motion of the field 
lines, and so the entire field-line stretching has to be attributed solely to the bipolar poloidal 
electric field. This illustrates the usefulness of this bipolar electric field as an important 
marker for two-fluid effects in the reconnection process. We also note that this electric field 
plays an important role in ion dynamics inside the reconnection layer. Namely, it is this 
field that is responsible for accelerating ions towards the midplane, leading to two counter- 
streaming ion beams and thus to an effective ion heating. The quadrupole toroidal magnetic 
field also plays an important role in ion dynamics as it diverts the two beams out of the 
layer (via the Lorentz force), thereby creating the expected stagnation-point pattern for the 
ion flow. 
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APPENDIX A: ALTERNATIVE DERIVATION OF THE ELECTRON VELOC- 
ITY FIELD IN THE SIMPLE X-POINT GEOMETRY 

In this appendix we present an alternative derivation of the poloidal electron velocity 
(and hence the toroidal magnetic field) for our simple X-point magnetic structure described 
by equations (J2J) and (|5jl. We make the same two basic assumptions: frozen- in law for the 
poloidal electron flow (which is not altered by the electron pressure) and incompressibility. 
Let us split the poloidal electron velocity field into two parts: parallel and perpendicular 
with respect to the poloidal magnetic field. According to equation ()3L)|). the perpendicular 
velocity is given by 

E z z x Bp i 

Vpoi,± = c — 2 . (Al) 

-Dpol 

In terms of the scaled variables introduced in Sec. Ill CI we can write the x and y compo- 
nents of this velocity as 



pol 
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v ±,y 



v ±,y 



\E\ 



R 2 

n pol 



im 1 - 



(A3) 



The parallel part of the poloidal velocity can be written as 



VpQi,|| = X(x, y) B pol . 



(A4) 



The total poloidal velocity field has to satisfy the incompressibility constraint, 



V • v = V • v 



pol, 



Vv 



pol, 



0. 



(A5) 



Using expressions (|A2|) — (|A3J) . we can write the divergence of the perpendicular poloidal 
velocity, to lowest order in S/L, as 



V ■ V po l,± = V • v 



,...u - d v v ±,y = -=r 



(A6) 



The divergence of the parallel velocity can be written as 



V ■ Vpoi,|| = V • (ABp i) = Bpoi • VA = Bpoi d h A , 



(A7) 



where L is the path length along a poloidal field line. Since the divergence of the total 
velocity must be zero, we get an expression for A: 



i„(s) 



A(*,x) 



(V- v 



pol,. 



B 



pol 



L f dx 

= -R- ( V -Vpou) 5- 
* -DO J £>x 



Using our expression (jA6|) for V ■ Vp^j,, we get 



Using equation (JJJ) for the field line shape, y(x, ty), we have 



L 

\&,x) = - — \E\ J 



dx 



j (2^ + x 2 ) 3 / 2 



L \E\ x 



(AS 



(A9) 



(A10) 
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Correspondingly, the components of the parallel velocity are 



m = HkL = = -Six, (All) 

^ - y-^-iy (A12) 

By combining this result with the components of the perpendicular velocity, we finally 

get 

\E\ _ 

v x = v\\ iX + v ±>x ~ v\\ >x = - x , (A13) 

\E\ ( \ \E\ 

v y = d iy + ^y = -Y\w + 1 ) = ~w Vj (A14) 

— in complete agreement with our calculation in Sec. Ill CI 



APPENDIX B: 2D STATIONARY IDEAL ELECTRON MHD: GENERAL FOR- 
MALISM 

In this section we describe a general formalism for analyzing a translationally-symmetric 
electron-MHD system in a steady state. We assume here that the electron density is uniform. 
In addition, we neglect any ion currents, both toroidal and poloidal; this assumption is valid 
if, for example, the ion temperature is negligible. We also neglect electron inertia; however, 
we do include an isotropic electron pressure in our equations. 

In full generality, the system is described by three vector fields: B, v* 6 ), and E, and a 
scalar electron pressure p e ; it is thus quite complicated. However, these fields are not all 
independent of each other. It turns out that the eMHD framework is so constraining that, 
with the help from the time-stationarity and translational symmetry conditions, the magnetic 
and electron velocity fields can be expressed in terms of only a single one-dimensional function 
and a constant. In the following, we outline how this is done. 
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First, from V ■ B = and the translational symmetry with respect to z, the magnetic 
field can be represented by two functions: the poloidal flux function "$?(x, y) (which in this 
section is measured from the X-point x = = y) and the toroidal field B z (x,y): 



B = BJ + B pol = B z z + V x [^z] = B z z + [W x z] . (Bl) 

The poloidal magnetic field components are 

B x = d v V; (B2) 

B y = -9 X *. (B3) 



Next, the electron velocity v* 6 ) is completely determined in terms of the magnetic field 
by Ampere's law: 

v (e) = -A = _ivxB, (B4) 
n e e D 

where D = Aim e e/c. In particular, 

W W = - 1 [V x [V x (*z)]], = 1 V 2 M/ , (B5) 
vgi = - \ [V x {B z z)\ = - 1 [VS, x 5] . (B6) 

Since the density is uniform, D is constant. 

Now let us turn to the electric field E. Generally speaking, the overall global magnetic 
configuration evolves as a result of reconnection. In particular, there is a continuous transfer 
of poloidal magnetic flux through the X-point and hence there is a non-zero toroidal electric 
field at that point. This electric field is a measure of the reconnection rate; it is inductive in 
nature. However, if the reconnection process is changing quasi-statically on the dynamical 
(Alfven) time scale, then locally, inside and around the layer, the magnetic field is essentially 
stationary, <9 t B = 0. Faraday's law then gives V x E = 0. Because of the translational 
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symmetry the toroidal electric field then has to be uniform; then, the poloidal electric field 
has to be potential: 

E z (x,y) = const; (B7) 
E P oi(z,2/) = -V(h(x,y). (B8) 

Thus, the full three-dimensional electric field is described by a constant inductive reconnec- 
tion field E z and a 2D electrostatic potential <j>2(x,y). 

This electric field is tied to the magnetic and velocity fields by the generalized Ohm law, 
i.e., the electron equation of motion. Neglecting the inertial and resistive terms but taking 
into account the Hall term and the electron pressure gradient term, this law becomes 

E + - [v^ xB] + ^ = 0. (B9) 

c n e e 

By taking the vector product of this equation with B, we get an expression for the perpen- 
dicular electron velocity: 

F, v R 

+ — • (BIO) 



(e) _ ExB c Vp e x B 

r± ~° n e e 5 2 



The first term in this equation is the ExB drift and the second term is the diamagnetic 
drift. Thus, the main effect of the pressure gradient term is the diamagnetic current. This 
current is in addition to that associated with the guiding center motion due to the ExB 
drift and needs to be included in the total current that one substitutes in Ampere's law. 

Provided that the electron density is uniform and the electron pressure tensor is isotropic, 
the diamagnetic currents do not lead to any substantial change in the mathematical structure 
of our formalism and are easily incorporated into our analysis. Indeed, combine the electric 
and pressure terms in Ohm's law into one by defining the modified electric field vector: 

E = E + —Vp e . (Bll) 

n e e 
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As we mentioned earlier, the poloidal electric field is potential, E po i = — V02- Then, since 
we also assume that the electron density is uniform in space, E po i is also potential, i.e., 

E = E z z- V0 2 = E z z- v(<h- — ), (B12) 

V n e ej 

and so vf } = c [E x B]/B 2 . 

This means that we can take into account the diamagnetic currents resulting from the 
electron pressure gradient simply by working with <p 2 and E instead of 02 and E. Note 
that one thus cannot really distinguish between the electrostatic potential and the electron 
pressure, and hence between the E x B drift and the diamagnetic drift, within the eMHD 
framework. This degeneracy however is not important; in particular, the magnetic and 
velocity fields can still be uniquely determined. 

The generalized Ohm law can now be written as 

E = - - [v (e > x B] . (B13) 
c 

Using equations ()B5j) - (jB6j) for and equation (jBljl for B, we can then express the com- 
ponents of E as 

E z = E z = — [VB Z x W] 2 = const , (B14) 
Epoi = " ^ (V 2 ^W + B Z VB Z ) . (B15) 

But expressing E po i in terms of the 2D potential 4>2(x, y) by equation (jB12|) . we see that 
(V 2 \E')V\1/ = Vg(x,y), where g = cDfa — B z /2. By taking the curl of this equation, we 
obtain V x (V 2 ^W) = V(V 2 ^) x W = 0, and hence B po i • V(V 2 *) = 0. That is, the 
toroidal current, j z ~ V 2 \l/, has to be constant along the field lines, and so must be a function 
of \l/ only: 

V 2 * = F(». (B16) 
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This equation is a necessary condition for the system to have a stationary solution, provided 
that the ion currents are neglected. This result is consistent with our earlier (Sec. IIII Aj) 
finding based on the analysis of the toroidal component of the eMHD magnetic induction 
equation, B pol • Vv z = v pol • VB Z . Notice that, in the main part of the paper we considered 
a general poloidal field configuration which could be supported by both ion and electron 
toroidal currents and thus equation ()B16|) did not need to be satisfied. However, the poloidal 
field structure of our specific example of Sec. Ill CI does satisfy this equation, which means 
that it could be produced purely by electron toroidal currents with no ion contribution. 

Once condition (|B16|) is satisfied, we can use the above formalism to compute, one by 
one, all the other electro-magnetic quantities. To do this for a given function F(^), one 
first solves the Poisson equation (|B16j) to find ^/(x,y), and then computes the poloidal 
magnetic field B pol and the toroidal velocity v z using equations (jBl|) and (jB5|) . Next, one 
uses equation (jB14|) . supplemented in the upstream region by the boundary condition B z (x = 
0,y) = [and by B z (x,y = 0) = in the downstream region], to calculate the toroidal 
magnetic field B z . Indeed, the meaning of this equation is that the rate of change of the 
toroidal field along a poloidal field line is equal to cDE z / B po \. 

Bpoi • VB Z = [W x z) ■ VB Z = [VB Z x W] ■ z = cDE z = const . (B17) 

Therefore, using the symmetry boundary condition at the y-axis, the toroidal field can be 
immediately obtained by integration along the field line: 

f dl i r dx 1 

*.(,, ,) = cDE z V(x, •) . cDE, j -g-*-^ = cDE , / ^-^ . (Big) 

Finally, one uses equation (jB6|) to determine the poloidal electron velocity and equation (jB15|) 
to determine the modified poloidal electric field. 
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This concludes the solution of the problem for a given poloidal field. The question of 
what determines the poloidal magnetic field structure lies beyond the scope of this paper. 
Here we just would like to remark that this structure is going to be affected by the pressure 
associated with the toroidal magnetic field. Specifically, since the toroidal field is strongest 
near the separatrix, it will push the flux surfaces apart, resulting in a weaker poloidal field 
near the separatrix. This effect has to be taken into account self-consistently. 
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FIG. 1: The basic idea of out-of-plane field generation. 
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FIG. 2: Simple X-point configuration. 
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FIG. 3: Bipolar poloidal electric field. 
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